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Limit  Cycle  of  the  Planar  Non-linear  Motion  / 

of  a  Finned  Projectile 

( 

Han  Zipeng 

[Abstract]  In  this  paper, "^the  planar  non-linear  motion  of  a 
finned  projectile  with  small  damp  and  weak  non-linear  static 
moment  has  been  studied.  The  stable  limit  cycle  is  used  to 
explain  the  self-oscillation  generated  by  some  projectiles.  The 
unstable  limit  cycle  is  used  to  investigate  the  causes  for  the 
fortuitous  fall  of  some  finned  projectiles.  In  order  to 
alleviate  the  effect  of  self -oscillation  and  to  eliminate  the 
fortuitous  fall  of  projectiles,  the  non-linear  portion  of  the 
damp  moment  and  that  of  lift  should  be  increased  and  the  non¬ 
linear  part  of  the  drag  should  be  decreased.-. 

1.  First  Approximation  Solution  of  the  Non-linear  Equation 
of  Motion  of  a  Finned  Projectile 

When  the  non-linear  aerodynamic  characteristics  are 
considered,  the  angular  equation  of  motion  of  a  finned  projectile 
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where  6-angle  of  attack;  v-velocity  of  the  center  of  mass  of  the 
projectile;  e-angle  of  Inclination  of  the  trajectory;  p-density 
of  air;  S-characteristic  area;  d-diameter  of  the  projectile;  A- 
equatorial  rotational  Inertia;  l-characteristic  length;  c  and 

Xq 

C  -constant  and  quadratic  coefficients  of  the  drag  coefficient; 

C  and  C  -linear  and  cubic  coefficients  of  the  lift  coefficient; 

Yo  ya 

m  and  m  -linear  and  cubic  coefficients  of  the  static  moment 

Zo  Zj 

coefficient;  and  m'  and  m'  -constant  and  quadratic  coeffici- 
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ent  in  the  equatorial  moment  coefficient. 

Based  on  various  aerodynamic  shapes  of  the  projectiles,  the 
position  of  the  center  of  gravity  and  the  Mach  number,  the 
algebraic  signs  of  b#  and  bj  may  have  the  following  four 
combinations^^^ : 

(1)  bo<0,  b,>0;  (2)  b,>0,  b,<0;  (3)  bo>0,  b,>0;  (4)  bo<0, 

ba  <0 

In  the  following,  we  used  the  van  der  pol  method^^^^^^^^^  to 
study  the  motion  of  a  finned  projectile  with  small  damp  (|bo/ 

•^k  |<<1)  and  weak  non-linear  static  moment  (  |k  6  [  <<  |k  1). 

The  majority  of  the  projectiles  belong  to  this  case. 

Let  us  leave  the  conservative  moment  k  6  in  equation  (1.1) 
on  the  left  sid^  and  move  the  remaining  moment  terms  to  the 
right.  Equtaion  (1.1)  can  be  written  in  the  following  form: 
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where 
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where  "j"  represents  taking  the  derivative  with  respect  to  the 
variable  t=(ijs  and  01=  \/-k  .  If  we  neglect  the  interference  term 

f(6,d)  of  small  damp  and  weak  non-linear  static  moment,  then  a 
linear  vibration  <j+6  =0  can  be  obtained.  Its  solution  is  a 
simple  harmonic  motion 

6  =  R  cos(t  +  ?)  (1.7) 

d=-R  sin(T  +  ♦)  (1.8) 


where 
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The  amplitude  R  and  the  phase  angle  f  are  determined  by  the 
initial  conditions  60  and  d*  (  which  is  the  derivative  of  arc 
length  s  with  respect  to  6).  When  the  effect  of  f(6,6)  is 
further  considered,  it  is  believed  that  the  solution  of  equation 
(1.1)  still  possesses  the  form  of  that  in  equations  (1.7)  and 
(1.8).  In  this  case,  however,  R  and  v  are  no  longer  constants. 
Instead,  they  are  functions  of  t.  a  and  6  can  be  derived  from 
equation  (1.7).  Then,  equations  (1.8)  and  (1.1)  are  used  to 
obtain  the  two  equations  relating  to  dR/dx  and  d^/dr.  Based  on 
these  equations,  we  can  solve 

+»).  “^^i»(t  ■«-V))*>(  T  +  ^  )  (1.10) 


Because  R  and  »  vary  for  slower  than  the  variable  C=t  + 
therefore,  in  the  first  approximation  the  right  side  of  the  two 
above  equations  can  be  considered  as  periodic  functions  of  the 
variable  C.  The  period  is  2ll.  Thus,  the  right  side  of  both 
equations  above  can  be  expanded  into  Fourier  series.  In 
addition,  the  effect  of  various  oscillation  terms  in  the  Fourier 
series  can  be  omitted  in  the  first  approximation.  Only  constant 
terms  remain.  It  is  then  substituted  into  the  expression  for 
f(6,6)  (i.e.,  equation  (1.6)).  The  following  result  is  obtained 
by  calculation 
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(1.13) 

2.  Phase  Plane  Analysis 

In  the  xoy  phase  plane,  the  position  of  the  phase  point  B  is 
determined  by  x=6  and  y= 6.  The  motion  of  the  phase  point  in  the 
phase  plane  reflects  the  oscillation  of  the  projectile.  Based  on 
equations  (1.7)  and  (1.8),  however,  the  position  of  the  phase 
point  B  can  also  be  determined  by  the  phase  radius  R  and  phase 
angle  t  +  9  (t)  (See  Figure  2.1). 

Based  on  van  der  pol's  method,  a  moving  phase  plane  oCn  is 
chosen,  which  rotates  clockwise  at  an  angular  velocity  w.  In 
this  moving  phase  plane,  the  polar  coordinates  of  the  phase  point 
B  happen  to  be  R  (t)  and  »(t).  Therefore,  equations  (1.12)  and 


(1.13)  describe  the  pattern  of  the  motion  of  the  phase  point  in 
the  moving  phase  plane. 

Let  us  first  discuss  the  case  of  linear  static  moment  where 
k  =0,  ♦(R)=0  and  •  Hence,  singular  points  are  located  along 
a  line  at  an  azimuth  9.  The  phase  radius  of  the  singular  point 
is  determined  by  F(R)=0. 


Figure  2.1  Dynamic  Phase  Plane  and  Fixed  Phase  Plane 

(1)  When  b, <0 ,  bj >0 ,  we  know  from  Equations  (1.12) 

(1.13)  that  there  are  two  singular  points  in  the  dynamic 
plane;  i.e.,  (0,0)  and  (Rp=2  V^bTTbJ ,9o ) •  When  0<R<Rp, 

the  radius  of  the  phase  point  continues  to  increase.  When  R>Rp , 

• 

R<0  and  the  phase  radius  continues  to  decrease.  Only  at  R=Rp, 

• 

R=0  and  the  phase  radius  remain  unchanged.  Therefore,  the 
second  singularity  is  a  stable  nodal  point  in  the  dynamic  phase 
plane.  Correspondingly,  this  singular  point  draws  a  circle  with 
a  radius  R^  in  the  fixed  phase  plane.  When  the  phase  point 
approaches  this  singular  point  radially  along  9o  in  the  dynamic 
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phase  plane,  a  spiral  can  be  drawn  in  the  fixed  phase  plane. 
Hence,  the  circle  is  the  stable  limit  cycle  in  the  fixed  phase 
plane.  In  this  case,  the  origin  is  an  unstable  focus.  (See 
Figure  2.2(a)). 


Figure  2.2  A  Stable  Limit  Cycle 


Figure  2.5  An  Unstable  Focus 


A  limit  cycle  is  an  isolated  closed  trajectory.  A  stable 
limit  cycle  represents  the  self-oscillation  of  the  projectile. 
Because  there  is  only  one  limit  cycle,  the  projectile  is  always 
excited  to  approach  the  same  steady  vibration  regardless  of  the 
initial  perturbation.  [See  Figure  2.2  (b),  (c)] 

6  =  RpCosCiDS  +  9o  )  (2.1) 

In  this  case,  the  motion  of  the  projectile  has  a  stable 
trajectory.  Because  k  =0  and  ^=0,  the  frequency  of  the  steady 
vibration  is  identical  to  that  of  the  linear  oscillator.  The 
formation  of  the  steady  vibration  and  the  magnitude  of  the 
amplitude  can  also  be  derived  based  on  equation  (1.12)  by  direct 
integration  and  by  letting  s-*®.  (The  derivation  is  omitted.) 

Mathematically,  when  k  =0,  b# <0  and  b* >0 ,  equation  (1.1) 
happens  to  be  the  van  der  pol  equation.  It  had  long  been  proven 
nathenaticallv  that  the  van  der  Pol  equation  has  a  stable  limit 
cycle.  This  limit  cycle,  however,  it  not  a  circle.  The  precise 


limit  cycle  with  a  medium  damp  (2bo /k  a;  1.0)  is  shown  as  curve  L 

Zo 


in  Figure  . 


Figure  2.6  Precise  Limit  Cycle  L  and  Approximate  Limit  Cycle 

L  (in  unit  of  y-\ 
a 

From  the  energy  point  of  view  we  can  further  comprehend  the 
self-oscillation  of  the  projectile  as  represented  by  the  limit 
cycle.  Let  us  multiply  both  sides  of  equation  (1.1)  by  6',  we 
get 

- 2  (6.+M*)(&'  )•  ( 2,.  2 ) 

2 

where  T  =  1/2(6')  represents  the  angular  kinetic  energy  of  the 

2  4 

projectile  and  V  =  -1/2  [k  6  +  (k  ]  represents  the 

Zq  Zj 

potential  energy  of  the  static  moment. 


T  +  V  represents  the  total  energy  of  the  angular  motion.  Based 
on  equation  (2.2)  we  can  see  that  the  angular  energy  increases 


Z  2 

when  6  <-bo /ba .  The  total  energy  decreases  when  6  >-bo /hz  .  If 
the  amount  of  energy  absorbed  by  the  projectile  in  a  cycle 
happens  to  replenish  the  energy  consumed,  then  the  projectile  can 
maintain  a  steady  vibration.  The  self-oscillation  represented  by 
the  motion  of  the  phase  point  along  the  limit  cycle  has  such 
characteristics . 

2  2  2 

Because  the  square  of  the  phase  radius,  R  =6  +(6'/\Ak  )  , 

is  proportional  to  the  total  energy  when  k  =0,  therefore,  the 

^2 

limit  cycle  radius  in  the  region  1 6 1  <'\/bo  /ba  (shaded  area  in  the 

figure)  gradually  increases  on  the  precise  limit  cycle  Lp  shown 

in  Figure  2.6.  Outside  this  region,  the  radius  is  gradually 

decreasing  because  the  energy  lost  in  the  segments  DE  and  FC  is 

compensated  by  that  gained  in  segments  EF  and  CD.  Therefore,  the 

projectile  can  stay  in  steady  oscillation. 

Self-oscillation  occurs  when  b# <0  and  b* >0  because  bo <0 

means  that  the  damp  is  negative  at  small  attack  angles.  In  this 

case  there  is  energy  input  in  the  angular  motion.  The  damp  is 

positive  at  large  angles  of  attack  (bjC  > | bo  |)  and  it  consumes 

energy  in  the  angular  motion.  The  energy  input  in  the  angular 

motion  is  obviously  converted  from  the  translational  kinetic 

energy  of  the  projectile.  This  can  be  qualitatively  explained  by 

2 

the  expression  of  bo  .  Because  bo  contains  a  term  (b  +gsin0/v  )<0 

(rising  arc),  it  acts  as  a  negative  damp  which  causes  the  angular 

energy  to  increase.  Based  on  the  equation  of  motion  of  the 

center  of  mass ,  (dv/d^-b^v  -gsine,  we  know  that  that  term  also 

makes  the  velocity  of  the  center  of  mass  and  the  translational 

2 

kinetic  energy  1/2mv  decrease.  Thus,  the  amount  of  energy 
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gained  angularly  Is  converted  from  the  translational  kinetic 
energy.  In  terms  of  the  physical  concept,  the  total  energy  of 
the  projectile  only  includes  the  angular  and  translational 
kinetic  energy.  Thus,  an  increase  in  the  angular  kinetic  energy 
can  only  be  converted  from  the  translational  kinetic  energy.  For 
projectiles  of  certain  aerodynamic  profiles,  the  drag  is 
particularly  large  because  the  equatorial  moment  and  lift  are 
small.  In  this  case,  the  situation  bg <0  may  occur. 

Consequently,  it  self-oscillates .  This  is  a  special  motion 
different  from  any  linear  motion. 

We  notice  that  mechanical  energy  conservation  exists  on  the  /66 
approximate  limit  cycle  (the  circle  R=Rp).  Consequently,  it  is 
not  possible  to  display  the  variation  of  energy  in  the  angular 
motion.  This  is  obviously  an  error  caused  by  the  first 
approximation.  From  Figure  2.6  we  can  also  see  that  the 
approximate  limit  cycle  is  more  or  less  the  average  position  of 
the  precise  limit  cycle.  The  smaller  b#  and  bg  become,  the 
closer  the  precise  limit  cycle  approaches  a  circle.  Thus,  the 
smaller  the  effect  of  this  error  becomes. 

(2)  When  bg >0  and  bg <0 ,  based  on  a  smaller  analysis  we  know 
that  the  R=Rp  circle  in  the  fixed  phase  plane  is  an  unstable 
limit  cycle.  The  origin  is  a  stable  focus.  Therefore,  when  the 
initial  condition  is  0<R  <Rp»  R<0 ,  the  phase  point  is  damped  to 
the  origin  [see  Figure  2.3  (a)].  The  corresponding  oscillation 
curve  is  shown  in  Figure  2.3(b).  The  motion  of  the  projectile  is 

4 

Stabilizing.  When  R>Rp i  R>0.  The  phase  point  is  far  away  from 


the  limit  cycle.  The  corresponding  oscillation  curve  is  shown  in 
Figure  2.3  (c).  The  motion  of  the  projectile  is  unstable. 

(3)  When  bo >0  and  b2<0,  only  the  origin  is  a  stable  nodal 
point.  In  addition,  we  always  have  R<0.  Hence,  the  motion  of 
the  projectile  is  stable  (See  Figure  2.4). 

(4)  When  bo <0,  b2<0,  we  always  have  R>0.  Hence,  the  motion 
of  the  projectile  is  unstable. 

Finally,  let  us  discuss  the  case  that  ^0.  From  equation 
(1.13)  we  know  that  The  phase  point  must  also  rotate  around 

the  origin  in  the  dynamic  phase  plane.  The  variation  of  the 
phase  radius,  however,  is  still  equation  (1.12).  Therefore,  the 
circle  Rp  =  2  V“t)o  /ha  is  still  the  limit  cycle.  Other  phase 
trajectories  are  elongated  clockwise  (k  <0)  or  squeezed 
counterclockwise  (k_  >0).  Correspondingly,  the  vibration  fre- 
quency  of  the  projectile  will  be  higher  and  lower  than  the 
frequency  of  the  linear  oscillator.  Furthermore,  it  causes  the 
frequency  to  be  related  to  the  amplitude.  However,  the  stability 
of  the  motion  is  not  affected. 

3.  Conclusions  and  Analysis 

Based  on  the  above  analysis  we  know  that  the  motion  of  a 
stable  static  finned  projectile  with  weak  non-linear  static 
moment  and  small  damp  is  becoming  stable  when  (1)  b,  >0,  ba>0. 

(2)  When  bo<0,  b2<0,  the  motion  of  the  projectile  is  unstable. 

(3)  When  b# <0  and  b2>0,  there  is  a  stable  limit  cycle  (the 
circle  R=Rp)  in  the  phase  plane.  In  this  case,  it  is  possible  to 
create  a  self-oscillation  at  an  amplitude  Rp.  (4)  When  b# >0  and 


ba<0,  there  is  an  unstable  limit  cycle  (the  circle  R=Rp)  in  the 
phase  plane.  If  R>Rp •  then  it  may  be  in  unstable  motion. 

In  the  design  experiment  of  finned  projectiles,  the  second 
case  described  must  undoubtedly  be  eliminated.  The  linear 
theory,  however,  neglects  the  non-linear  damp  term  ba .  It 
believes  that  as  long  as  bo>0,  a  stable  static  finned  projectile 
must  be  in  stable  motion.  Most  of  the  finned  projectiles  were 
also  thus  designed.  In  reality,  however,  bj  does  exist.  The 
following  two  situations  may  appear.  One  is  that  ba >0  over  the 
entire  trajectory  so  that  non-linear  stability  coincides  with 
linear  stability.  Thus,  the  effect  of  bj  is  masked.  The  other 
is  ba <0  over  the  entire  trajectory  or  on  a  portion  of  the 
trajectory  to  produce  an  unstable  limit  cycle  in  the  phase  plane. 
When  the  initial  disturbance  is  excessive  or  when  encountering  a 
gust  in  flight,  the  phase  point  may  jump  out  of  the  unstable 
limit  cycle,  i.e.,  R=/6  +(6'/w)  >Rp'  Hence,  the  motion  becomes 
unstable,  leading  to  the  fortuitous  fall  of  projectiles.  Because 
the  probability  of  having  an  extraordinarily  large  initial 
disturbance  or  gust  is  very  small,  especially  on  trajectories 
where  bj <0  occurs  only  in  certain  sections,  the  fall  of 
projectile  due  to  the  unstable  limit  cycle  is  accidental. 

As  for  fallen  finned  projectiles,  we  should  inspect  their  b, 
values  at  various  Mach  numbers.  If  bj<0,  we  should  alter  the 
aerodynamic  profile  and  layout  to  make  b,  >0  or  to  minimize  |  bjJ  . 
Because  the  radius  of  the  unstable  limit  cycle  gets  larger  with 
smaller  jb*  | ,  the  range  of  initial  conditions  which  satisfies 


.■V  m” 
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that  of  a  stable  motion  is  also  getting  wider.  When  IbjI-O, 

Rp  Then,  any  initial  condition  cannot  cause  unstable  motion. 

In  addition,  certain  odd-shaped  projectiles  with  small 
equatorial  moment  and  lift  and  large  drag  may  encounter  the 
situation  that  bo <0  and  b* >0  at  some  Mach  numbers.  In  this  case 
the  projectile  will  remain  in  self-oscillation  at  certain  attack 
angles.  If  the  amplitude  of  this  oscillation  is  too  large,  the 
projectile  is  flying  at  a  large  angle  of  attack  over  a  long 
period  of  time,  leading  to  decreasing  range  and  deteriorating 
degree  of  concentration.  Therefore,  the  aerodynamic  profiles  of 
this  type  of  projectiles  should  also  be  changed.  The  value  bj  /67 
should  be  increased  to  reduce  the  radius  of  the  stable  limit 
cycle  in  order  to  allow  the  projectile  to  fly  at  a  small  angle  of 
attack. 
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Akstract 


In  this  ptptT,  non-linear  motion  of  finned  projectile  with  small  damp 
and  weak  non-linear  static  moment  has  been  stadied.  Using  the  stable  limit 
cycle*  it  has  been  explained  that  some  finned  projectiles  probably  produce 
self-oscillation. 

Using  the  anstable  limit  cycle*  tke  causes  of  the  fortuitous  fall  of  some 
finned  projectiles  have  been  discussed.  In  order  to  decrease  the  influence 
of  the  self-oscilation  and  eliminate  the  fortuitous  fall  of  the  finned  pro¬ 
jectile*  the  non-linear  part  of  the  damping  moment  and  lift  must  be  inc¬ 
reased*  and  the  non-linear  part  of  the  drag  must  be  decreased. 


